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ABSTRACT 


With the huge advance in artificial intelligence and the rapid development of 
intelligent swarm algorithms, the exploration of facility layout problem 
(FLP) with its non-deterministic polynomial-time (NP-Hard) nature has 
gained much more attention. The squirrel search algorithm is one of the 
swarm algorithms that is known for its effective gliding feature that provides 
cheap exploration of lengthy distances. In this work, Msqrl algorithm is 
presented as a modification of squirrel search algorithm to be capable of 
handling permutation-specific FLP. The modification is done by introducing 
two new operators: Msqrl-Exchange and Msqrl-Winter. It is used to 
investigate the effectiveness in finding acceptable solutions to variable-size, 
single-row FLPs in a fast and efficient manner. Tests included small and 
large benchmark instances for comparisons. Outcomes show that Msqrl was 
able to improved quite a few previously found solutions by acting efficiently 


and converging rapidly to solutions. It outperformed both semidefinite 
programming and cuckoo optimization in finding optimal solutions in an 
acceptable number of iterations and relatively small population sizes. 
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1. INTRODUCTION 

Facility location or facility layout planning problems (FLP) outline a vital subdivision of 
optimization problems, which are responsible of optimally positioning facilities to minimalize an objective 
function based either on effort, cost, time, or distance. Facilities are elements, machines or workstations that 
usually enable activity functioning, and their layouts resemble configurations essential for performing certain 
operations. In the framework of production systems, the main purpose is minimizing overall costs for 
material handling involved in the movement of materials among sections. It has been found that nearly up to 
20-50% of the entire operational manufacturing cost goes to materials handling. FLP has many practical 
application areas such as planning layouts for gardens, stores, hospitals, offices, and plants [1]. 

In view of the high demand in manufacturing industries for reduced costs, efficient facility layout 
designs are sought. This problem is a recognized combinatorial optimization problem, with its 
non-deterministic polynomial-time (NP-Hard) complexity, it is considered to be difficult to solve [2]. Finding 
the exact solution for a certain FLP problem means enumerating all possible arrangements of layouts, with n 
facilities the number of permutations is n!. That is why only small problems can be treated using exact 
algorithms. Such exact methods usually require excessive memory are considered be expensive in terms 
of computations [2], [3]. That is way, other intelligent methods are being investigated such as metaheuristics 
and swarm intelligence. These methods can develop fast approximated near-optimal solutions for larger 
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problems. One of the recently developed swarm algorithms is the squirrel search algorithm, it replicates the 
behavior of southern flying squirrels with their economical way of movement (locomotion) that gives them 
the ability to modify lift and drag. As squirrels fly in the forest trying to find food sources, they glide 
changing locations on trees in search of nuts to eat (acorn) and nuts for winter storage (hickory). Researchers 
have several objectives to think about when dealing with FLPs, such as reducing the cost of material 
handling, decreasing material flow, dropping down the overall distance traveled for the material, and even 
refining the entire adjacency among facilities. 

Through the last decades, authors have addressed various types of this problem, as in 2005, 
Anjos et al. constructed a relaxed semidefinite programming (SDP) method to find lower bounds on the 
optimal values by using the natural symmetry of the problem [2]. In 2008, Khilwani et al. used 
Psycho-Clonal algorithm in relation to the dynamic system characteristics and operational constraints 
to solve the FLP, it was found that the suggested algorithm outperformed others [4]. In 2009, Amaral 
suggested a new lower bound for the single-row FLP through the optimization of a linear program over the 
partial description at hand with cutting planes [5]. In 2010, Ohmori et al. proposed solving FLP using swarm 
intelligence employing particle swarm optimization (PSO) [6]. In 2011, Datta et al. used genetic algorithms 
(GA) that is based on permutations to solve single-row FLPs [7]. Prasad et al. in 2014, designed manufacturing 
plant layouts using computerized relative allocation of facilities (CRAFT) and developed a JAVA program 
to design the optimum layout using STEP file as input [8]. In 2015, Matai used a modified simulated annealing 
to solve the multi objective facility layout problem [9]. Goncalves and Resende presented an encoding 
of unequal area FLPs using GA grounded on biased randomized number [10]. In 2016, Anjos and Vieira 
initiated a framework to solve unequal area layouts using a combination of two models based on mathematical 
optimization [11]. Park and Seo suggested a two-step heuristic algorithm in 2017, the first step uses 
a construction method which sequentially places facilities, and the second step enhances that of the first step. 
Results indicated that the proposed algorithm yielded equivalent quality of layouts within much shorter time 
when compared with former work [12]. Zhao et al. incorporated GA with Levy Flight for improving mutations, 
they also used an island model genetic algorithm with varying approaches of decoding every island [13]. 
In 2018, Kalita and Datta studied the constrained single-row facility layout problem using GA with large 
problem instances [14]. Moslemipour et al. used a hybridization among ant colony, clonal selection and 
simulated annealing for large FLPs [15]. Feng and Chae, made the effort of solving fixed-shape rectangular 
layouts through the use of integer linear programming, they employed new constraints to maximize flow of 
material between the closed facilities [16]. In 2019, Lin and Yingjie presented survival signature for layout 
problems along with optimizing solutions using GA [17]. Also in 2019, Kim and Chae employed the 
monarch butterfly optimization with a slicing tree to denote the layout [18]. Kalita and Datta in 2020, 
attempted to solve the constrained Single-row FLP using GA with repair procedures, the method was applied 
to a set of different-size FLP instances [19]. Tongur et al. employed migrating bird optimization, tabu search 
and simulated annealing to find layouts for a big scaled university hospital [20]. Moreover in 2020, Meskar 
and Eshghi a used a lower bound to solve the generalized single-row FLP using cuckoo optimization 
algorithm [21]. Lately in 2021, Ingole and Singh used biogeography-based optimization to find solutions for 
unequal area (fixed/flexible) FLPs [22]. 

In this work, the squirrel search algorithm is to be modified for tackling the problem of varied 
length, single-row facility layout problems. The modified algorithm (Msqrl) is investigated using benchmark 
flow instances less than 30 [7], and others between 60 and 80 [2], two operators are introduced to modify the 
algorithm (Msqrl-Exchange and Msqrl-Winter). Results are compared with pervious work done using the 
same data instances to show the effectiveness of the modifications carried out in this work. 


2. RESEARCH METHOD 
2.1. Problem formulation for facility layout planning 
Facility layout planning represents a very interesting real-world class of problems, a variety of 
observational studies can be found in the literature. Going through the literature, the definition of FLP is 
given as “a common industrial problem in which the objective is to configure facilities, so as to minimize the 
cost of transporting materials between them” [23]. As for optimization problems, definite optimal solutions 
are sought, so there will be a need for sufficient methods capable of recognizing optimal solutions [24]. 
Layouts for facilities can be stated using an arrangement of entities necessary to carry out certain 
responsibilities such as producing profits or services delivery. Each facility represents a unit object for the 
implementation and execution of any task. This ranges from machine shops to departments or warehouses [23]. 
For optimization problems, definite optimal solutions are sought. So, there will be a need for some 
sufficient methods capable of recognizing at least one optimal solution [24]. Mathematically, the objective is 
minimizing the distance-based measure given in (1). 
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where flow; is the flow quantity between department i and j, dj; is the direct straight distance between them, 
ci is the cost of moving one unit from i to J. 

Given two dissimilar facilities f; and fj, and a permutation p. The value of Dp(i, j) is the summation 
of the lengths for facilities from the center of f; to the center of f; (taking into consideration 1f the facilities are 
placed on the right or left side of the other facility). Solving the problem requires finding a permutation that 
gives the minimized sum of distances between all facilities pairs. Mathematically [2]: 


M Zi low ‘Dy Ci j) (2) 


In this work, a special case of facility layout planning is studied for comparison reasons. It is called 
the linear single-row facility layout problem. This problem involves placing facilities of various lengths 
on a straight line. 


2.2. Squirrel search algorithm 

Squirrels fly among trees for foraging and fulfilling their needs. They glide from one tree to another 
seeking resources of food in autumn, during this time, they vary their location to discover diverse areas of the 
forest. When weather temperatures get sufficiently hot, squirrels can find their everyday energy requirements 
fast enough by instantly consuming acorns whenever found. Having completed their daily supply, they begin 
a new quest for a food source that is best for winter, that is hickory nuts. Storing nuts can aid in sustaining 
energy needs in severe meteorological conditions and cuts the expenses of foraging journeys. In cold seasons, 
forest trees drop their leaf cover exposing squirrels to predators, so they tend to be less energetic with no 
hibernation. Such a process keeps repeating itself during the life of a squirrel, establishing the basis of the 
algorithm. Squirrel search algorithm presumes (n) flying squirrels, each is found on a separate tree in a forest 
and separately looks for nuts using accessible food resources optimally with some behavior of dynamic 
foraging. In the forest, three types of trees exist: normal, oak (acorns) and hickory trees, assuming that the 
forest includes three oak and only one hickory tree [25]. 


2.2.1. The basic algorithm 
Starting with a randomly originated vector S containing locations for N squirrels, they begin gliding 
in the search space to modify their locations. Where Si j denotes iteration ¢ with the i squirrel in the j” 


dimension. Next, is the S matrix showing the representation of squirrels: 


SN, SN,2 _ Sw,d 


A fitness function f (S i j) is used to evaluate the locations of squirrels, the squirrel with the minimum 
fitness value is said to be on the hickory tree (storing food) and denoted as (SÉ jh» the subsequent three best 
fitness values are for squirrels presumed to be on acorn trees (consuming normal food) and are denoted 
as (Si ja The rest of flying squirrels are on normal trees (Si j)n with no resource of food. As they change 
locations, a number of the flying squirrels move in the direction of the hickory tree having achieved their 
daily needs, others move to acorn trees to get their daily requirements. 

The mechanism of gliding is defined by the lift L computed in (3) and drag D force calculated in (4). 
Flying squirrels glide at steady speed and drop at angle (@g) horizontally with lift-to-drag ratio. 
The descending angle of gliding is shown in (5) [25]. 


L = 1 / 2pCLV2S (3) 


where p is the air density and equals 1.204 kgm”. Cy is the lift coefficient usually set in the range 
0.675 < C_< 1.5, Vis the speed and equals 5.25 mst. S is the body surface area and equals 154 cm’. 


D = 1 / 2 pV2SCD (4) 


where Cp is the frictional drag coefficient, usually set to 0.6. 
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D 
p = arctan(7) (5) 


by substitution, we get: 


1/2 pV?SCp 


= arctan 
P ( 1/2pC,V2S 


) (6) 
p = arctan (P) (7) 


Meaning that the only Cp and Cz are controlling (@) and are both constants. Now, the constant of 
gliding distance G4 can be computed as in (8). 


Ga = Hioss /(SCF * tan (@)) (8) 


where Hioss is the height loss after gliding set to 8 m, ScF is a scaling factor set between 16 and 37. 

Flying squirrels update their locations based upon three factors. These are their previous location, 
gliding distance, and probability (Ppa) that a predator is close. They move according to three circumstances: [25] 
— Flying squirrels on acorn trees can head to hickory tree to preserve optimal food resources by (9). 


if rand 2 Ppa (SijJa + Ga XC x (Sin — Sia) 
Otherwise Random location 


(Si) =| (9) 


where C is a constant of gliding to governor the traversing of the search space, rand is a randomly generated 

number in [0, 1]. 

— Squirrels residing on normal trees are able to glide in the direction of acorn nut trees to satisfy daily 
consumptions. The resulting location of squirrels can be updated as in (10). 


if rand > Ppa (Sij)n + Ga X C X ((SÉ;)a — (SÉ;)n) 


(10) 
Otherwise Random location 


(Sn =| 


— For squirrels that has consumed acorn nuts and are now on normal trees, they will possibly fly gliding 
for the hickory tree to start storing food for winter. Hence, updated locations can be obtained from (11). 


if rand > Pya(Sijn + Ga X C x (Sin — (Sin) 
Otherwise Random location 


(SD =| (11) 


2.2.2. Observing seasons 

In the cold winter season, flying squirrels lose substantial amount of heat. Owing to their raised 
temperature and small body size, searching for food becomes costly and also risky because of predators. 
Unlike autumn, winter influences and reduces their activity. Due to its importance, a procedure for observing 
and controlling the season check is add to avoid getting stuck in local optima. The season constant Se 
is calculate using (12). 


7 : 2 
Sc = JEk-1(Sax = Sny) (12) 


where i is the number of squirrels on acorn trees, d is the dimension length of a squirrel. 
Next, seasonal check is carried out (SE < Smin), here Smin is the minimum value of a season constant 
calculated as [25]: 


10E76 
Smin — 365t/(tm/2.5) (13) 


where f¢ is the current iteration, tm is the maximum iteration. 

Higher values of Smin encourage exploration, whereas lesser values improve the ability of 
exploitation for an algorithm. There should always be a suitable stability between them. If checking for 
seasons becomes true, meaning winter has ended, flying squirrels that were unable to discover optimal food 
for winter are repositioned in a random fashion. Figure 1 gives the detailed flowchart of the squirrel search 
algorithm and the process of finding optimal solutions. 
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2.3. Modification of squirrel search 

Due to the nature of the FLP, a variation is carried out to modify the squirrel searching 
methodology. This modification addresses the main movement steps of squirrels along with other 
adaptations. The process is mainly grounded on discovering good-enough solutions through some certain 
tactics to renew locations of squirrels, thus producing new populations targeting for global optimums. Each 
squirrel represents a sequence vector of N facilities numbered form 1 to N, every single one of these 
sequences is a suggested layout solution. 


Set no. of Squirrels=N, max iteration= MaxlIter, current 
iteration: Iter=1, initiate first population randomly 





Calculation of fitness values for each individual Population 


locations of squirrels are sorted ascending on fitness values 


Y 






Season Check satisfied? 
Yes 


Random repositioning of flying squirrels after winter 
Update (Smin) the minimum of the 
lter=lter+1 season constant using Eq. (13) 
Yes 
No 
Stopping Criterion Iter=MaxIter? 


Exit with the final optimal location of the squirrel on the hickory tree 





Figure 1. Flowchart for squirrel search algorithm 


The modifications suggested in this work involve three issues. These are: correctly creating initial solutions, 
updating locations of squirrels, and monitoring season check. They are conducted as follows: 


2.3.1. Initial solution 

Initializing the individuals of the first population of flying squirrels is performed randomly. Here the 
permutation characteristic of FLP should be kept in mind. Accordingly, the first population is randomly 
generated from the sequence of facilities ranging from | to number of facilities N. 


2.3.2. Locations update for squirrels 

In the original algorithm, squirrels update their location in three ways: 

— Squirrels on acorn trees head for the hickory nut tree, 
— Squirrels on normal trees head for direction of acorn trees, and 
— Squirrels on normal trees head for the hickory nut tree. 

This is done using (9), (10), and (11), these equations yield incorrect solutions when applied to FLP. 
In order to preserve the correctness of the produced solutions a new operation Msqrl-Exchange (©) 
is introduced to replace subtraction and gliding. This operation uses a mechanism to find a set of location 
pairs if exchanged, new acceptable sequences will be obtained as solutions to FLP. 

The Msgqrl-Exchange operation takes two squirrel sequences as input and performs a subtraction like 
process to imitate the movement of squirrel from one type of trees to the other. This is done by finding pairs 
of locations to exchange in the second sequence without repetition of locations (locations that are already 
been chosen for exchange are neglected). The process is explained with a fully detailed example in Figure 2. 
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As for the two constants Gd and C used in the original work. Both constants have no impact on 


updating locations of squirrels using the Msqrl-Exchange operation. They can be omitted without affecting the 
procedure of the algorithm. 


Locations 1 2 3 4 5 6 7 8 9 
suieta = | 3[ 4] 5} 1] 9] 8] 6] 7 | 2 | 
= > << 
w 
Se 


Squirrel 2 


All pairs of location: (1,9) (2,4) (3,1) (4.8) (5.5) (6,6) (7.3) (8,7) (9.2) 


Neglectedas 
locations 2,9 have 


Neglected as location I 


AY 
has already been chosen No need to 


exchange 
E already been chosen 


Neglected as location 4 Neglected as location 8 has 
has already been chosen already been chosen 


| 


Pairs of locations to exchange in Squirrel 2 are: (1,9) (2,4) (7,3) 


(3.7) 
(2.4) (1.9) 


Process of exchange: 
Locations 


1 2 3 4 5 6 7 8 9 

suirel2 | § | 2 | 6 {4 | ots | 7a] | 
4 

swied?  (aļaj|7|2jo|s|s|i]s] 


Figure 2. An illustration of the proposed Msqrl-Exchange process with an example 





In achieving all modifications mentioned so far, (9), (10), and (11), become: 


(S a = : rand = Ppa (Sipa + (Spx © (Si j)a) (14) 
4 Otherwise Random location 

(Sa = A rand > Ppa Gn + (CSE; Ja © (Si j)n) (15) 
n Otherwise Random location 

Cc, = {i rand = Ppa Gai + (Sin © (Si j)n) (16) 
@ Otherwise Random location 


2.3.3. Adjustment for winter season 

The season constant Sc measures the closeness of squirrels on acorn trees to the hickory tree (best 
achieved solution in the population). It should be possible to accomplish the same effect using another new 
adjustment operator called Msqrl-Winter (H) instead of subtraction. So, having two squirrels with n facilities 
each, only locations with the same value are counted, this is illustrated as follows: 


Squirrel 1: 5 2 6 4 7 8 1 9 J 
H 
Squirrel 2: 5 2 6 7 8 9 1 4 3 
1 1 1 0 0 0 1 0 Lej 


That is, squirrels 1 and 2 have five similar locations out of nine. Hence, (12) can be rewritten as in (17). 
This way, we have a correct measure for the season constant ready for use with FLP. 
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d A 
ac Irak HSh,k) % 


Se = X% i C1 (17) 


where ac is the no. of acorn trees, d is the dimension, C/ is a constant set to 0.01. 


3. RESULTS AND DISCUSSION 

This section involves testing the effectiveness of the modified squirrel search algorithm Msaqrl 
in solving FLPs. Flow matrix instances with lengths of facilities are taken from the literature, small problem 
instances (<30) [7] and larger numbers between (60 and 80) [2] are considered for comparisons. The distance 
matrix is calculated for each generated layout in the course of the methodology. Solutions are compared with 
semidefinite programming (SDP) [2], permutation-based GA (PGA) [7], and cuckoo optimization algorithm 
(COA) [21]. 


3.1. Parameter setting 

For small numbers of facilities, setting the parameter involves taking (50) squirrels as the population 
size taken in the original work [25]. For food resources, (1) tree is for hickory and (3) trees are for acorn nuts, 
and the rest (46) trees are normal trees without nuts. As for larger facility numbers, (100, 150, and 200) 
squirrels are assigned: (1) hickory tree, (6, 9, and 12) acorn trees, and the remaining are normal trees. The 
probability of a predator appearing Ppa is (0.1). The number of facilities, maximum iteration and the total 
number of possibilities n! are shown in Table 1. 


Table 1. Parameters for tests along with the number of possibilities 
No. of facilities | No. of squirrels __ No. of iterations Possibilities n/ 


Small problems 5 50 5 120 
8 50 5 40,320 

10 50 10 3,628,800 

11 50 10 39,916,800 

20 50 30 2.432902e+18 

30 50 100 2.6525286e+32 

Large problems 60 100 500 8.3209871e+81 
70 150 1000 1.197857e+100 

75 150 2500 2.480914e+109 

80 200 4000 7.156946e+118 


3.2. Tests and comparisons 

In order to inspect the proficiency of the modified algorithm, two tests are carried out. The first one 
uses six traditional instances with identified optimal solutions, these problems have no clearance between the 
facilities [2]. The obtained results in Table 2 indicate the capability of Msqrl algorithm in achieving the same 
results with a small number of iterations and squirrels. 


Table 2. Comparison of the Msqrl using facility numbers < 30 (Test1) 


rae Methods 
No. of facilities SDP [2] PGA [7] Msarl 


5 151.0 151.0 151.0 
8 2324.5 2324.5 2324.5 
10 2781.5 2781.5 2781.5 
11 6933.5 6933.5 6933.5 
20 15549.0 15549.0 15549.0 
30 44965.0 44965.0 44965.0 


On the other hand, the second test uses 20 instances for larger problems having no optimal solution. 
These instances include 60, 70, 75 and 80 facilities, and are divided into 4 groups of 5 instances each. These 
were originally created in random by [2] and are frequently used for comparisons in the literature. Table 3 
illustrates the results gained by former methods and the modified algorithm Msqrl. As can be seen from the 
results, the Msqrl algorithm was capable of enhancing 5 solutions out of 20 (shown in bold along with the 
layout sequence of the enhanced solutions). As for the remaining 15 instances, results are the same with those 
gained by PGA which are the best results among other. 

Results of Testl together with Test2 signify the outstanding performance of Msqrl. The algorithm 
was successful in surpassing both SDP and COA and finding better solutions. This was achieved 
using a relatively small number of iterations and population sizes of squirrels. 
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Table 3. Comparison of Msqrl algorithm using facility numbers 60-80 (Test2) 


P SDP [2] COA [21] PGA [7] Msqrl Best layout (for improved solutions only) 


60 1,493,704 1,477,834 1,477,834 1,477,834 

843,644 841,814 841,792 841,790 53-10-13-41-29-50-30-9-15-42-8-14-34-32-40-33-28-21-56-37-25- 
5-24-1-54-39-17-46-45-18-36-52-43-6-59-51-26-55-58-27-16-60- 
38-47-7-49-3 | -3-12-48-23-2-4-19-22-57-35-44-11-20 

656,272.5 648,651.5 648,337.5 648,337 

405,433 398,481 398,468 398,424 13-20-3 1-19-55-57-10-39-7-44-24-22-36-33-56-11-21-6-23-17-46- 
40-59-53-5-3-45-50-2-35-14-27-15-18-58-8-16-60-29-26-42-37-54- 
47-30-5 1-32-48-52-9-49-38-28-1-4-41-34-25-43-12 

319,501 318,855 318,805 318,805 


70 1,543,098 1,528,760 = 1,528,621 1,528,621 

1,494,182 1,441,504 1,441,028 1,441,028 

1,524,171.5 1,519,578.5 1,518,993.5  1,518,993.5 

974,856 968,796 968,796 968,796 

4,230,912.5 4,218,017.5 4,218,017.5 4,218,003.5 67-28-17-44-64-32-69-33-18-13-52-47-25-57-48-45-55-15-59-42- 
66-3-46-60-34-29-7-1 1-53-21-4-30-36-5-12-70-20-35-14-24-56-43- 
41-61-26-8-49-27-63-65-16-54-6-10-39-22-50-37-58-40-2-51-31- 
62-23-68-1-9-38-19 


75 2,399,583.5 23,958,845 2,393,456.5 2,393,456.5 
4,348,544 4,321,381 4,321,190 4,321,190 
1,295,085 1,248,664 1,248,537 1,248,537 
3,949,276.5 3,942,749.5 3,941,891.5 3,941,816.5 21-17-9-49-72-38-69-68-27-25-6-51-61-74-11-31-4-71-43-54-59- 
58-23-26-66-63-56-48-55-19-2-67-1-32-13-45-24-18-34-73-44-53- 
28-52-16-65-35-64-3-12-46-39-33-40-29-41-20-57-22-47-30-70-8- 
37-62-42-10-75-7-50-15-14-5-60-36 
1,816,455 1,792,038 1,791,408 1,791,408 


80 2,138,083.5 2,069,192.5 2,069,097.5 2,069,097.5 
1,939,938 1,921,202 1,921,177 1;921,177 
3,332,421 3,251,435 3,251,368 3,251,368 
3,773,429 3,747,548 3,746,515 3,746,515 
1,611,495 1,588,963 1,588,901 1,588,892 25-17-48-41-58-60-69-78-13-49-42-72-12-29-16-5 1-14-64-26-28- 
24-15-56-5-57-50-46-23-63-40-45-20-80-4-18-77-39-32-3-59-65- 
55-66-7 1-35-53-43-6-44-54-75-79-34-30-62-33-36-76-9-19-3 1-8- 
70-1-22-11-74-37-10-68-67-21-73-61-27-38-52-47-7-2 


4. CONCLUSION 

In this work, the squirrel search algorithm was successfully modified with a newly introduced 
operators (Msqrl-Exchange and Msqrl-Winter) to be applied correctly for solving single-row facility layout 
problems. Small and large facility numbers were considered in the testing process and comparisons were 
made with instances taken from the literature to show the efficiency of the modified algorithm Msgqrl. In view 
of all the tests and results achieved in this work, the modified algorithm has a fast convergence to solutions, 
it was able to find better solutions to 5 of the used instances and it and outperformed both SDP and COA. 
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